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Abstract 



Using the quark mass density- and temperature-dependent model, we have 
studied the properties of the dibaryon systems. The binding energy, radius 
and mean lifetime of (f2f2) +and are given. We find the dibaryons 

(f2$7) +, are metastable at zero temperature, but the strong 

decay channel for (f2f2) + opens when temperature arrives at 129.3 MeV. It 
is shown that our results are in good agreement with those given by the chiral 
5(3) quark model. 

PACS number: 12.39.Ki, 14.20. Jn, 11.10. Wx, 24.85.+p 
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I. INTRODUCTION 



Since Jaffe predicted the existence of H particle [1], the study of dibaryon systems has 
attracted more and more attentions [2-14]. A lot of theoretical models, for example, the 
constituent quark model in which the effective degrees of freedom are constituent quarks 
and gluons [12]; the quark-mesons exchange model in which the effective degrees of freedom 
are constituent quarks and Goldstone bosons [9,13]; the chiral SU(3) quark model [6-8,14], 
... had been employed to predict the bound states of dibaryon systems. Many particles, 
namely, H dihyperon [1-3], d! particle [5], (fifi) + and ITS [6-8,14], AA [9,11] have been 
suggested to be candidates for being observed in the experiment. Unfortunately, the con- 
vincing experimental evidence for the existence of these particles has not yet been found. 
Then it is essential to confirm which dibaryon would be the best one to be detected in the 
experiment. Therefore, it is necessary to address the above dibaryons from different point of 
views, different models and different treatments, because it can open the field of vision and 
choose the best detectable candidate after comparing the results given by different models. 
This is the first motivation to employ our quark mass density- and temperature-dependent 
(QMDTD) model [15-18] to study the dibaryon system. 

The second motivation is to extend the study of dihyperons to finite temperature. The 
strangeness enhancement was observed in relativistic heavy ion collision (RHIC) recently 
[19,20]. The dihyperon can be found from the extreme condition of RHIC easily than 
that from the usual one. According to the thermal model of RHIC, a fireball with high 
temperature and high density is formed. Employing a suitable model to study the stability 
and the thermodynamic properties of dibaryon at finite temperature can help us to learn 
the detectable possibility of dihyperon. 

In previous papers, by means of the QMDTD model, we studied the thermodynamic 
behavior and stability of A = 5 and A = 10 strangelets at finite temperature where A 
is the baryon number [15,18]. We found that the stable strangelets only exist in the high 
strangeness module and high negative charge region. For A = 2 dibaryon system, (f2f2) + has 
the maximum strangeness module and maximum negative charge. It is of interest to study 
A = 2 dibaryon system by using QMDTD model. 

The organization of this paper is as follows. In next section, we give the main formulas 
of the QMDTD model. The results including stability, binding energy, radius, lifetime for 
dibaryon systems at zero temperature are presented in Sec. III. In Sec. IV, we will discuss 
the stability of dibaryon via possible strong and weak decays at finite temperature. The last 
section contains a summary. 

II. QMDTD MODEL 

QMDTD model is a non-permanent quark confinement model [16] because it is based on 
the Friedberg-Lee model [21]. It has been used to study the thermodynamic properties of 
strangelets [15,17,18] and the photo strange star [22]. The detail of this model can be found 
in [15-18]. Here we only write down the main steps which is necessary for calculating the 
thermodynamic quantities. 

According to the QMDTD model, the masses of u, d quarks and strange quarks (and the 
corresponding anti-quarks) are given by [15-18] 
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where n# is the baryon number density, itl s q is the current mass of the strange quark and 
B(T) is the vacuum energy density which satisfies 



B(T) = B 



T 



TV 



, < T < T c , 



B(T) = 0, T > T c , 



(3) 
(4) 



where -Bo is the vacuum energy density inside the bag (bag constant) at zero temperature , 
T c =170 MeV is the critical temperature of quark deconfinement phase transition, and a, 
b are two adjust parameters. The values of a, b are fixed as 



a = 0.65, b= -0.35. 



(5) 



Using the same argument as that of Ref. [6], in RHIC experiments, the fireball expands 
until freeze-out. The production of a large amount of particles and violent collisions between 
them in the fireball finally lead to the thermal and chemical equilibrium, soon after it reaches 
high temperature and high density. The thermodynamic potential in the fireball reads [6] 
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ln(l + 
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where % stands for u, d, s (or u, d, s ) quarks, fa is the corresponding chemical potential 
(for antiparticle fa = —fa). Si{k) = \Jmf + k 2 is the single particle energy and is mass 

for quarks and antiquarks. — ^ - is the density of states for various flavor quarks. For a 

dibaryon, its geometry can be taken as a sphere. The density of states of a spherical cavity 
has been calculated in our previous paper, it reads [23]: 



Ni(k) = Ai(kR) 3 + Bi(kR) 2 + Q(kR), 



(7) 



A,= 



2^ 
9tt' 



(8) 



B, 



9i_ 
2tt 



1 + 



mi 
k 



tan 1 



~k 



7T 
2 



(9) 



^ Oj f 1 ( k mA , k nk 

Ci = ¥-\-+\— + -^) tan" 1 

\rrii 



m 
- \ ~~r~~ 



1.45 



9i 



2;r |3 \ ///, k I mi 2rrii J \kj 3.42 (™± _ 6.5) +100 

where ^ is the total degeneracy. 
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Since the quark mass rrii depends on the density, the thermodynamic potential £1 is not 
only a function of temperature, volume and chemical potential, but also of density. The 
expression of total energy density reads [15-18] 



t an 
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where % = A/V is the baryon number density and A is the baryon number of the strangelet, 
4 

V = g 71 "-^ 3 1S the volume of the strangelet. The number density of each particle can be 
obtained by means of 



1 dtt 
V dm 



(12) 
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At finite temperature, we must include the contributions of the anti-particles, therefore, the 
baryon number for % quark is given by 



AiV; = (m - ni) x V 



l 
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The strangeness number \S\ of the strangelet reads 

\S\ = AN S: 

the baryon number A of the strangelet satisfies 

A = ^(AN u + AN d + AN s ). 
The electric charge Z of the strangelet is 



(13) 



(14) 



(15) 



Z = ^AN U - ^AN d - ^AN S . (16) 
At finite temperature, the stability condition of strangelets for the bag radius R reads 

where the free energy F of strangelet is 

F = E-TS, (18) 

E = eV is the total energy, and 



dT 



(19) 



is the entropy. For the dibaryon system, baryon number A, the strangeness number \S\, 
and electric charge Z are given. Solving above equations self-consistently, we obtain the 
thermodynamic quantities of dibaryon systems. 
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III. RESULT AT ZERO TEMPERATURE 



A. Stability 

Firstly, we use the same arguments as that of Ref. [24,15] to study the stability 
of dibaryons at zero temperature. For a dibaryon Q(A,\S\,Z), its baryon number A, 
strangeness number \S\ and electric charge Z are conserved in the strong process. A general 
expression of a two-body strong baryon decay for a dibaryon bound state can be written as 

Q(2, \S\,Z) - x(l, \S X \,Z X ) + y(l, \S\ - \S X \,Z - Z x ), (20) 

where x and y stand for baryons. For example, in the case of (f2f2) + , x and y are free Q, 
respectively. At zero temperature, the process (20) is allowed if the energy balance of the 
corresponding reaction satisfies 

E(2,\S\,Z)>m x + m y , (21) 

where E stands for the energy of the dibaryon bound state. 

In the weak process, the baryon number A and the electric charge Z are conserved, but 
the strangeness number \S\ is not conserved. For the weak decay, = ±1. Therefore, a 
general expression of a two-body weak decay for dibaryon can be written as 

Q(2, \S\,Z) - x(l, \S X \,Z X ) + y(l, \S\ - \S X \ ± 1, Z - Z x ), . (22) 

This process is allowed at zero temperature if the energy balance of the corresponding 
reaction satisfies 

E(2,\S\,Z) > m x + m y , (23) 

For example, the two-body weak decay to pure hadrons channel for (f2f2) + is 

(tttt) 0+ ^tt + Z-. (24) 

In addition, (f2fi) + can also decay to another dibaryon and a hadron via weak processes 

(nn) 0+ -> (QE°) +7T" (25) 

(nn) 0+ -> (HE-) +7T°, (26) 

provided 

E(2,6,-2)>E(2,5,-l)+m^, (27) 

£(2,6,-2) > E(2,5,-2)+m n o, (28) 

respectively. On the other hand, the three-body weak decay channels for (f2f2) + particle 
read 

(nn) + ^n + A + K~, (29) 

(nn) 0+ ^n + E° + 7r-, (30) 

(nn) 0+ ^n + E- + 7r°. (31) 
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These channels can open provided the energy E(2, \S\, Z) is large than the sum of masses 
of the corresponding three particles. 
We choose the following area 

N s - N, > 0, (32) 
Z > -A, (33) 
\S\+Z<2A. (34) 

to investigate the stability of dibaryon systems. By using the definition of Schaffner-Bielich 
et. al [24], a dibaryon is called unstable if the inequality (21) is satisfied because in this 
case the strong decay occurs. A dibaryon is called metastable if the inequality (21) is not 
satisfied but at least one of the weak decay inequalities is satisfied, in this case the weak 
decay occurs. Other cases in which all inequalities be broken are called absolutely stable. 

Noting that the free energy F reduces to the energy E at zero temperature, we can 
study the stability of dibaryon system by using the formulas of Sec. II for QMDTD model. 
Our result is shown in Fig. 1 where open circles stand for the unstable dibaryons and filled 
circles for the metastable dibaryons. We find no dibaryon can be absolutely stable at zero 
temperature. But three dibaryons, (fifi) + , (fiS - ) and (f2H°) corresponding to Q(2,6, —2), 
<5(2,5,— 2) and (5(2,5,-1), respectively, are metastable. For (f2f2) + particle, its strong 
decay channel and a weak channel (29) are forbidden, and the other weak channels, (24)- 
(25), (30) and (31), are permitted. The above conclusions are the same as that of Ref. [6-8], 
though results of [6-8] are established on the chiral SU(3) quark model and ours on the 
QMDTD model. 



B. Binding energy and stable radius 



Secondly, we study the binding energy and the stable radius for above metastable 
dibaryons. The binding energy of (f2f2) + is 



B(nn) + = 2xm n - E(nn) + . 

The stable radius R of (f2f2) + is determined by the condition 

5E(nn) 0+ 



5R 



= 



(35) 



(36) 



at zero temperature. Obviously, the binding energy E(QQ) + depend on the values of the 
two parameters Bq and m S Q. The permitted ranges of Bq and m S Q are called "stability 
window" of the strange quark matter. We have calculated this stability window in Ref. 
[16] for QMDTD model. The binding energy and stable radius of (fif2) + state for fixed 
m s0 = 150 MeV and different B are shown in Table 1: 



B Q MeV fnT a 


m s0 MeV 


B(nn) 0+ MeV 


i?fm 


170 


150 


162.7 


1.13 


175 


150 


143.7 


1.12 


180 


150 


125.2 


1.11 


182.5 


150 


116.1 


1.11 


190 


150 


89.2 


1.10 



(Table 1) 
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We find that the value of binding energy of (f2f2) + calculated by our model is in good 
agreement with that by chiral SU(3) quark model. Hereafter we choose parameters Bq = 
182.5 MeV fm -3 and m s0 = 150 MeV, because the binding energy 116.1 MeV is just equals 
to that predicted in Ref. [6]. With these parameters, we draw the energy of (f2f2) + state 
as a function of radius in Fig. 2. It is found the stable radius R — 1.11 fm is a little bigger 
than 0.84 fm given by chiral SU(3) quark model. 

To compare our result with chiral SU(3) quark model furthermore, we also calculate the 
binding energy of (f2S _ ) dibaryon 

B(QE~) =mn + m s - - E(nE~). (37) 

and its stable radius. The results are shown in Table 2: 



B MeV far 3 


m s0 MeV 


B(QEr) MeV 


R fm 


170 


150 


114.0 


1.13 


175 


150 


96.5 


1.12 


180 


150 


79.4 


1.12 


182.5 


150 


71.0 


1.11 


190 


150 


46.2 


1.10 



(Table 2) 



We find that the values of binding energies of the (f2S~) is also consistent with the chiral 
SU(3) quark model [7]. The difference between (f2f2) + and (fiS~) is about several tens 
MeV. In particular, we hope to emphasize that the QMDTD model gives the same conclusion 
as that of chiral SU(3) quark model: the dibaryon (f2f2) + is a deeply bound state, and its 
binding energy is maximum in dibaryon systems. 



C. Mean lifetime 

The key for the detectability of the dibaryon is to study its mean life time. Instead of 
the calculation of two-body decay Feynman diagram [6], we employ a decay formula which 
is first given by Chin and Kerman [25] to estimate the mean life time of (fif2) + and (fiS~) 
via weak decay. The decay formula is 

l/ r = [GVs/1927r 3 ] sin 2 9 c F(z), (38) 

with 

F(z) = 1 - 8z + 8z 3 - z 4 - I2z 2 In z, (39) 

z = /i 2 //! 2 , (40) 

where the Cabibbo angle is given by sin6> c ~ 0.22. With Eqs. (38)-(40), we can roughly 
estimate order-of-magnitude of the dibaryon's mean life time. The coupling constant G is 
determined as follows: using Eq. (38) to calculate the decay mean lifetime of fl for following 
channels: 

tt^A + K~, (41) 
n^Z° + ir-, (42) 
n -> 2° + tt° (43) 
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and fit the value with the Q mean lifetime 8.21 x 10~ n s to determine coupling constant 
G. Then use this coupling constant G to calculate the lifetime of the (f2f2) + and (QS~). 
With the parameters B = 182.5 MeV fm -3 , m s0 = 150 MeV, we obtain the mean lifetime 
is 4.7 x KT 11 s for (fifi) 0+ and 6.2 x KT 11 s for (QS~), respectively. 



IV. RESULT AT FINITE TEMPERATURE 

Now we turn to discuss the thermodynamic properties, especially, the stability of 
dibaryons at finite temperature. Since according to the chiral 577(3) quark model and 
QMDTD model, the best stable candidate of dibaryons is (fifi) + at zero temperature, we 
focus our attention only on (f2f2) + in this section. Instead of the energy at zero temper- 
ature, we calculate the free energy density first. Given the strangeness number, baryon 
number and electric charge, the free energy and stable radius for dibaryons and baryons can 
be obtained self-consistent ly by using the formulas in Sec. II. 

As was pointed out by [6-8] and sec. Ill, the strong decay channel 

(tttt) 0+ ^tt + tt, (44) 

is forbidden because the binding energy of (f2f2) + is large. But at high temperature, a 
metastable dibaryon can absorb enough energy from the hot environment to overcome its 
binding energy and open the strong decay channel. To study the possibility of strong decay, 
we treat the (Qfi) + and Q as a six quarks cluster with Q(2, 6, —2) and a three quarks cluster 
with (5(1,3, —1) respectively, and calculate their free energy density. Define 

Af = f(nn) 0+ - 2f(n) (45) 

as the difference between the free energy density of (f2f2) + and two separate Q system, the 
curve of Af vs. temperature T is shown in Fig. 3. We find Af increases with temperature. 
It becomes zero when temperature arrives at T = 129.3 MeV. When T > T , Af > 0, the 
strong decay channel opens. This result is reasonable if we notice that the temperature T 
is a little bigger than the binding energy of (f2f2) + , 116.1 MeV. 

Since the strong decay of (f2f2) + will open and (£7f2) + becomes unstable when T > 
T , the favourable detectability of (fif2) + must be taken in the regions < T < 129.3 
MeV. This temperature regions seem wide enough for detectability in RHIC experiments. 
Of course, in order to help the design of a detectable experiment, it is necessary to get 
more information about (f2fi) + in this metastable regions. For this purpose, we study the 
temperature dependence of the radius and the mean lifetime of (fifi) + in < T < T 
regions. 

The curves of stable radius R vs. T for (f2fi) + is shown in Fig. 4. We see that the 
radius increases from 1.11 fm to 1.89 fm when temperature increases from to 129.3 MeV. 
the volume of the bag expands when temperature increases. 

Finally, we still use Eqs. (38)- (40) to study the temperature effect on the mean lifetime. 
Instead of Fermi energies of quarks at zero temperature in these equations, we use the tem- 
perature dependent chemical potentials fi u (T) and /i s (T) to estimate the order of magnitude 
for the mean lifetime of (On) + . We find that the order of magnitude is not less than 10~ n 
s in the whole metastable temperature regions. 
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V. SUMMARY 



In summary, we have studied thermodynamic properties of dibaryons by using QMDTD 
model. At zero temperature, we find that only three dibaryons, namely, (f2f2) + , (fiS - ) and 
are metastable. At least one weak decay channel opens for these states and all the 
strong decay channels close. We calculate the binding energy, stable radius as well as the 
mean lifetime of (f2fi) + and (J1S~), it is found that our results are in good agreement with 
those given by chiral SU(3) quark model. Using two essentially different models and get 
the same conclusion, this result strongly compels us to support (f2f2) + is the most favorable 
candidate of dibaryon for detectability, because its binding energy is maximum. 

Extending our investigation to finite temperature, we find the strong decay channel will 
open and (f2f2)o+ will become unstable when temperature arrives at To = 129.3 MeV. We 
also calculate the radius and the lifetime of (f2f2) + at finite temperature. 
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VII. FIGURE CAPTIONS 

Figure 1. The electric charge Z as a function of the strangeness number \S\ for unstable 
strangelets (open circles), metastable strangelets (filled circles) with baryon number A = 2 
at zero temperature. 

Figure 2. The energy per baryon number E/A as a function of the radius R for 
(f2f2) + particle. 

Figure 3. The difference of free energy density between (f2f2) + particle and two free f2 
particles Af as a function of temperature T. 

Figure 4. The stable radius R as a function of the temperature T for (f2f2) + particle. 

VIII. TABLE CAPTIONS 

Table 1. The binding energy and stable radius of (f2f2) + particle. 
Table 2. The binding energy and stable radius of (f2H~) particle. 
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